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Research  on  the  GPS  Platform-Based 
ESG  North  Finder  Integration  System 


Chen  Jiabin 

(Department  of  Automatic  Control,  Nanjing  University 
of  Aeronautics  and  Astronautics,  Nanjing,  210016) 


Abstract;  The  GPS  platform-based  ESG  north  finder 
integration  scheme,  which  is  a  new  scheme  to  achieve 
north  seeking  on  its  own  and  is  of  practical  value,  is 
presented  in  this  paper.  It  has  the  advantages  of 
better  mobility,  high  accuracy,  quick  alignment, 
identification  of  the  ESG  drift  models  in  real  time. 

Key  words;  GPS,  platform  ESG  north  finder,  azimuth 
angle 


1 .  Introduction 


In  modern  land  combat  technology,  emphasis  has  been 
increasingly  placed  on  mobility  and  flexibility  because  raising 
the  mobility  of  armed  forces  and  impact  point  hit  rate  becomes 
one  of  the  key  measures  in  self  protection.  Owing  to  the  high 
cost  and  powerful  kill  ability  of  a  missile,  accurate  positioning 
and  direction  finding  is  expected  to  be  conducted  for  the  weapon 
system  before  a  missile  is  launched  so  as  to  accurately  destroy 
the  potential  targets.  At  present,  the  azimuth  alignment  adopted 
for  flying  missiles  uses  launch  pad  alignment  and  optical 
alignment,  while  ballistic  missiles  generally  employ  an  optical 
pointing  technique.  The  launch  pad  alignment,  however,  can 
provide  only  a  lower  precision  with  an  azimuth  alignment  error 
reaching  1-10^  while  optical  alignment,  though  reaching  a 
certain  precision  at  an  angular  second  level,  is  poor  in  mobility 
and  takes  as  long  as  0.5-lh  to  accomplish  an  alignment,  suitable 
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only  for  a  launch  site  with  a  known  geographic  location  and 
target  azimuth,  and  it  is  apt  to  expose  itself. 

To  raise  the  concealment  ability  and  mobility,  missiles  are 
expected  to  perform  north  finding  on  their  own  and  with  a  higher 
precision  at  a  launch  site  with  an  unknown  geographic  location  in 
as  short  a  time  as  possible.  A  platform  static  electricity  gyro 
north  finder  plan  presented  in  the  References  [2]  possesses  the 
capability  of  realizing  self  north  finding  in  a  short  time  but 
requires  the  geographic  location  of  the  launch  vehicle  to  be 
known.  To  solve  this  problem,  other  auxiliary  positioning 
equipment  has  to  be  applied,  such  as  GPS,  etc.  A  high  precision 
satellite  navigation  system  vigorously  developed  in  the  United 
States,  GPS  can  provide  various  users  with  information  of  high 
precision  three-dimensional  location,  three-dimensional  velocity 
and  accurate  time,  and  is  widely  applied  in  aviation,  aerospace, 
sea  voyage,  vehicle  navigation  and  geodesy.  Once  GPS  is  combined 
with  the  platform  static  electricity  gyro  north  finder  to  form  an 
integrated  system,  it  can  realize  accurate  positioning  and 
direction  finding  with  a  mobile  missile  launch  pad. 

This  paper  presents  a  brief  introduction  of  the  principle  of 
the  GPS/platform  static  electricity  gyro  north  finder,  which,  by 
using  the  approximate  non-linear  filtering  theory,  can  identify 
the  azimuth  angle  of  the  system,  the  location  of  the  launch  pad 
and  the  gyro  drift  model  coefficient.  Also,  it  can  carry  out  a 
theoretical  analysis  of  the  azimuth  angle  error  generated  by  the 
latitude  error.  The  study  shows  that  this  integration  system, 
with  a  higher  north  finding  precision,  serves  as  a  new  plan  of 
realizing  azimuth  alignment  during  a  mobile  missile  launch  and 
has  applications  in  engineering. 

2.  Principle  of  the  GPS/Platform  Static  Electricity  Gyro  North 
Finder  Integration  System 

The  principle  of  the  GPS/platform  static  electricity  gyro 
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north  finder  integration  system  is  shown  in  Fig.  1.  The  static 
electricity  gyroscope  is  mounted  on  a  servo  rotating  base  plate 
with  its  gyro  shell  rotating  around  the  rotor  momentum  matrix 
axis  at  an  angular  velocity  <i>^j  0  and  a  respectively  are  corners 
of  the  outer  ring  axis  and  inner  ring  axis  of  the  rotating  floor; 
0=a=O  at  the  initial  instant,  Ug  is  the  information  output  by  GPS 
receiver;  y  is  azimuth  angle.  When  the  alignment  time  is  no 
more  than  20min,  the  variation  law  of  and  a  caused  by  the 
Earth's  rotation  is  [2]: 


1^  fD  — oj,]  [ZT]  fsinyl  ^  fo-l 

y'u  ojy+'^wy  n, 

where  (Oj  and  02  north  and  vertical  components  of  the  Earth's 
rotation  angular  velocity  o,,  <ih=at.cos<p,  <a2=<a.siD'p,  a>,j 
where  Oy  is  gyro  drift.  From  the  References  [6],  the 
expressions  of  and  Wy  can  be  derived  as: 
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where  S=sin  Ott,  S=cos  Wbt,  A=[aia2a3a4a5bib2b3b4b5]'^  is  constant 
state  vector;  *«  and  u  are  Guassian  white  noise  of  zero  mean 
value. 
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Fig.  1.  Principle  of  GPS/Platform 
Static  Electricity  Gyro  North  Finder 

Key: 

(a)  Top  signal 

(b)  Equator  signal 

(c)  Control  circuit 

(d)  Rotating  floor  servo  circuit 

(e)  Corner  collection  circuit 

( f )  Computer 

1.  Static  electricity  gyroscope  shell 

2 .  Gyroscope  rotor 

3 .  Servo  electric  machine 

4.  Induction  synchronizer 

5.  Photoelectric  sensor 

6 .  Outer  ring  axis 

7 .  Inner  ring  axis 


Equations  (1)  and  (2)  indicate  that  the  variation  law  of  0 
and  a  is  associated  with  user's  latitude  (p,  azimuth  angle  y  and 
gyro  drift.  Also,  the  simultaneous  Eqs.  (1)  and  (2)  take  the 
state  of  the  GPS  receiver  into  consideration.  In  this  case,  the 
state  equation  of  the  integration  system  can  be  derived  as 
follows : 

x(t)f[X(t)]+[(t)W(t)]  (3) 

where  X(  t  )  — [0QtY/^1^2^3®'A^5^1^2^3^4^5^g]^/  ^  C^)  L*«  3  ! 
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are  the  state  vector  and  noise  of  GPS  receiver,  f[X(t)], 
[(t)]  are  given  different  values  for  different  integration  modes. 
The  observation  equation  of  the  system  is 

Y(t)=H[X(t)]+V(t)  (4) 

where  Y(t)  is  the  swept  observation  space,  V(t)  is  observation 
noise.  The  specific  forms  of  H[X(t)]  and  V(t)  will  be  detailed 
later. 

If  the  action  of  and  Oy  is  omitted,  the  azimuth  error  Sy 
caused  by  latitude  error  6(p  can  be  derived  from  Eqs.  (1)  and  (2) 
as 

Obviously,  in  domestic  regions,  the  selected  tp  value 
usually  will  not  be  too  small,  and  6(p  will  generate  a  larger  Sy 
when  y  is  larger  and  therefore,  it  is  necessary  to  conduct 
accurate  positioning  for  the  launch  pad. 

With  a  higher  positioning  precision,  GPS  is  an  ideal  system 
used  for  positioning.  With  a  horizontal  positioning  error  100m 
at  present,  this  system,  if  combined  with  the  platform  static 
electricity  gyro  north  finder,  will  be  extremely  useful  in 
increasing  north  finding  precision.  Different  GPS  receiver 
models  output  different  information.  When  the  output  information 
is  the  user's  longitude  and  latitude  and  altitude,  the 
implications  of  parameters  in  the  state  equation  (3)  and 
observation  equation  (4)  respectively  are; 


X(j,)  =  [_e  ay  ai  02  03  04  osbi  6?  6364  65  (pX 


/Cx(0]= 


r<i>2a— <i)iSiny +0  i5+o  jC + a(a  jC' -f-o  +  2^0  BiSC 

— <i>iCosy4-i>i54-i>2C+«(i5C*+6<5*)-f  2a  b,SC 
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(8) 

(9) 


—s  -c 

r(0=  C  -3 
.  Oi2xt  . 

>"(0  =  [K,F.Kj' 

where  Vo,  and  V,  are  Guassian  white  noise  of  zero  mean  value, 
while  other  parameters  have  the  same  implications  as  above. 

When  the  GPS  receiver  outputs  pseudo  range  information  and 
user  attitude  can  be  ignored,  the  implications  of  parameters  in 
state  equation  (3)  and  observation  equation  (4)  respectively  are 

X (j.')—[9ay  ai  aiO)  a,  as  b\  bib 
iuja— (i>ismv  4-a  i;S +a  aC 
q>20  ct)jCosy 
/LX(0]=  o,,^i 

I 

f  -s  -a  0  0  ■ 

c  —5  0  0 

r(t)=  ■ 

0  0  10 

.  0  0  0  1  J 

r.v  (i)  j= r^o  (>i  i>i"'p,y 

F(0  =  [P%  F,F,  V,—V,~j 


where  t^  and  t^^  respectively  are  the  range  error  when  the  user 
clock  error  is  equivalent  to  clock  frequency  drift  error,  and  the 
range  variation  rate;  A  is  user  longitude;  p(i=l,2. . . ,n)  is 
pseudo  range  between  user  and  satellite  i;  n  is  nvimber  of  visible 
quantities;  Vi(i=l,2. .  •  ,n)  is  pseudo  range  observation  noise; 
other  parameters  have  the  same  implications  as  above.  The 
expression  of  pseudo  range  p  is 


(10) 

(11) 
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A  —  ^  .cos  cos  A)  ^  -r  (y  ^ — fi  .cos  9»  s  inX )  “  +  (Z., — R,sin(p)*  7'* + /, 

(18) 

where  (X,i,  Yji  and  Zji)  is  the  location  of  the  satellite  i  in 
Earth's  coordinates,  which  can  be  calculated  from  the  sidereal 
data  provided  by  GPS  receiver. 


The  non-linear  time-varied  system  made  of  the  state  equation 
and  observation  equation  of  the  integration  system  has  no 
available  analytical  solutions  and  should  be  analyzed  by  using 
the  numerical  calculation  method.  In  the  study  of  the 
integration  system  with  the  approximately  non-linear  filtering 
theory,  the  initial  azimuth  angle  during  the  missile  launch,  the 
coefficients  of  gyro  drift  model  and  the  accurate  location 
information  of  the  missile  launch  pad  can  be  identified.  The 
parameter  identification  process  is  as  follows: 


Xt+i/i  —  Xijk~rf  (-X»/»)T «T ^  (X./i)/ (^Xm')T\/ 2 

Kt+x  —  P  i+iP  i+x/],HIj.x-\-Ri+i') 

x»+i/4+i=Xi+iA+^Lk+iCyi+i  H  (x»+i/»)] 
P»+l/»+l  =  [^~X»+lfl’i+l]P  »+!/» 


where .4 (A-.A)=afCx (<)]/& (OUO)=Xk A. 


aril) 


^  {.4[x(0]/[x(0])  •  ^  r,  =  riu),  X.+ia+i 


,=1' 


N"’ 


To  is  sampling  time;  recursion  initial  value  is  Xo/o/  Po/c 


(19) 

(20) 
(21) 
(22) 
(23) 


3 .  System  Simulation  and  Its  Result  Analysis 


Suppose  a  missile  launch  pad  arrives  at  an  unknown  launch 
site  with  true  longitude  and  latitude  A=118°,  (p=32°,  while  a 

target  under  attack  is  stationary,  located  at  north  to  east  20®. 
A  GPS  receiver  will  output  longitude  and  latitude  and  altitude 
info3nnation;  recursion  initial  value ; 
the  element  values  on  the  diagonal  line  are  respectively  listed 
in  Tables  1  and  2;  other  elements  are  zero. 


7 


Table  1 . 


Element  Values  on  P,  diagonal  line 


(oA*) 

(*)= 

i 

1  2  3 

4  5  6  7  8  9  10  n 

12  13  14 

0.5=  4= 

0.  001= 

10^ 

Key;  (1)  Unit 


Table  2.  Element 

Values 

on 

0fc,  R^.  Diagonal  Line 

(o/iO 

(.)2 

1 

1  2 

1 

2  3 

0.  001= 

0.  5* 

0.  3=  3.  5= 

Key:  (1)  Unit 


Table  3  lists  the  accurate  coefficient  values  of  the  gyro 
drift  model:  sampling  time  To=0.5S;  azimuth  alignment  time  is 
20min;  gyro  actual  body  rotation  angle  velocity  (Ob=Tr/h;  the 
forms  of  the  parameters  in  the  state  equation  and  observation 
equation  of  the  system  are  selected  from  Eq.  (6).  Based  on  the 
system  performance  analysis,  the  stable  estimated  values  of 
azimuth  angle  and  latitude  are  computed  as; 

y  =  20.  003°  9>=  32.0001° 

The  stable  estimated  values  of  the  coefficients  a^, 
bi(i=l,2,3,4,5)  of  the  gyro  drift  model  are  listed  in  Table  3. 
Figs.  2-4  show  the  variation  curves  in  the 

process  of  identification. 
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Table  3.  Coefficients  of  Gyro  Drift  Model 


o/h 

o/h 

i 

1 

2  3 

A 

5 

i 

1 

2 

3 

4 

5 

ai 

0.1 

0.2  0.4 

0.4 

a4 

C4 

0. 106 

0.193 

0.  325 

0.  306 

0.  510 

hi 

0.6 

0.7  0.5 

0,5 

0.5 

hi 

0.592 

0.  6^0 

0.  422 

0.  581 

0.  565 

Key:  (1)  Unit 
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The  system  simulation  result  shows  that  based  on  the 
location  information  output  by  the  GPS  receiver,  the  azimuth 
angle  and  the  location  of  the  missile  launch  vehicle  at  the 
initial  missile  pointing  can  be  more  accurately  estimated  by 
using  Kalman  filtering  theory  with  an  azimuth  error  around  10" 
and  latitude  error  around  0.4",  indicating  a  rather  high  north 
finding  precision.  In  addition,  this  integration  plan  also 
allows  a  real-time  identification  of  the  gyro  drift  model  to 
solve  the  problem  of  poor  repetition  of  the  drift  model 
coefficients  when  the  static  electricity  gyros  start  one  by  one. 
Undoubtedly,  this  is  extremely  helpful  in  raising  the  north 
finding  precision  of  the  system.  Therefore,  the  GPS/platform 
static  electricity  gyro  north  finder  plan  presented  in  this  paper 
appears  to  be  a  new  plan  of  realizing  self  north  finding  during 
mobile  missile  launch,  and  it  is  also  applicable  to  other 
occasions  which  need  accurate  north  finding. 
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